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Abstract We present a new version of a sample 
of galaxies from the Revised Flat Galaxy Catalogue 
(RFGC), which have redshift and Hi line width data. 
We also give the parameters of the collective motion 
model determined upon this sample. The considered 
models of motion include the dipole (bulk flow), the 
quadrupole (cosmic shear) and the octupole compo- 
nents. We also considered higher-order multipoles. In 
all cases the obtained parameters matched the ACDM 
cosmology. 

Keywords galaxies: kinematics and dynamics; galax- 
ies: distances and redshifts; galaxies: spiral; methods: 
statistical 



1 Introduction 

The distribution of matter density is inhomogeneous in 
the region of the Universe, limited to about 100h~^ Mpc,| 
which contains several superclusters and voids. A ga- 
laxy, besides the cosmological expansion, is also at- 
tracted to the regions with greater density. Due to this 
fact, the galaxies are involved in a large-scale collective 
motion on the background of Hubble expansion. Inves- 
tigation of such a motion is important since it allows 
plotting the distribution of matter in the surrounding 
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region of the Universe and comparing this distribution 
with the distribution of luminous matter. 

This is especially important due to the fact that the 
most serious challenges to the standard ACDM cos- 
mology are posed by the inconsistencies in the estima- 
tions of the velocity of the bulk motion. The ACDM 
model estimates it at the level ^ 250 kms~^ at the scale 
lOOfe-^Mpc. However, i n some studies, for e xample 
Lauer fc PostmanI (Il994l): IWatkins et all (l2009l): Feld- 



man et al. ( 20101 ) ■ the obtained values were larger by a 
factor o f 2. The largest valu e of 416 ± 79kms~^ was 
given bv lFeldman et all (|2010f ). Other results, including 
our own, give smaller values, which are consistent with 
the predictions of the ACDM cosmology. 

In addition to redshifts we require independent es- 
timations of distances to the galaxies. Since we deal 
only with spiral galaxies, we use the Tully-Fisher re- 
lation to determine these distances. In its common 
version it relates the intrinsic luminosity of a galaxy 
with its velocity width (the amplitude of its rotation 
curve). However, we use the 'Hi line width - linear di- 
ameter' variant of the Tully-Fisher relation, which does 
not require photometric data. It is different from the 
common version and the data are also processed in a 
different fashion. If these two variants give similar re- 
sults, this is a strong evidence of their correctness. For 
this reason, it is important to continue the research of 
the collective motions of galaxies, even if it can not 
boast the best accuracy or exceptionally large depth. 
Of course, different versions of the Tully-Fisher rela- 
tion can not be regarded as independent methods, but 
they can potentially give very different results. 

Due to large errors in determination of distances to 
galaxies, which are caused both by measurement errors 
and by the intrinsic uncertainty of the Tully-Fisher re- 
lation, it is necessary to compile large samples, and to 
pay special attention to the adequate processing of the 
data. Also, the choice of the model of the collective mo- 
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tion strongly affects the results. For low-depth samples 
it was possible to use the simple bulk motion model. 
However, this model is inadequate for deeper samples 
and more complex models should be used. These mod- 
els contain higher-order multipoles and take into ac- 
count additional effects. 

We study these motio ns using the Revised Flat Ga- 
laxy Catalogue (RFGC. [Xarachentsev et aL 1999 ). It 
contains data about N = 4236 galaxies including the 
information on the following parameters: Right As- 
cension and Declination for the epochs J2000.0 and 
B1950.0, galactic longitude and latitude, major and mi- 
nor blue and red diameters in arcminutes in the POSS- 
I diameter system, morphological type of the spiral 
galaxies according to the Hubble classification, index 
of the mean surface brightness (I - high, IV - very low) 
and some other parameters, which are not used in this 
article. More detailed description of the catalog ue can 
be found in the paper ([Karachentsev et al.lll999f ). 



2 Data used 

We used the sample containing 1720 RFGC galaxies. 
The data for 59 galaxies were rejected due to large de- 
viations from the Tully-Fisher relation. Thus, the re- 
sulting sample contained 1661 galaxies. In addition to 
RFGC data we used the data on the radial velocities 
of the galaxies and H i line width at the 50% level. As 
usual, all velocities were converted to the CMB frame 
and H i line widths were corrected for intrinsic absorp- 
tion and turbulence. The new data were taken from 

3 releases 



the SFI-I— I- II survey (ISpringob et al 
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of th e ALFALFA survey (GioyaneUi et al.ll2007l: Kent 



et al. 200 81: iMartin et al.ll2009f ) and the 40ALFA sur- 
vey ((Havnes et al.ll201lr ). which contains about 40% of 
data to enter the final version of the ALFALFA sur- 
vey. Some results of RFGC galaxies observation at the 
Effelsberg radio telescope were taken from the papers 
(jMitronova et al.ll2005l : IXudrva et"aDl2009l ). From ah 
other sources only one mea surement was incl uded in 
this sample from the paper (jKovac et al.l 120091 ) . which 
contains the results of a blind survey in the Canes Ve- 
natici region. Thus, the new sample is based on much 
more homogeneous data than the previous one, which 
is an additional benefit. 

As a result, we added the data about 42 new galax- 
ies, 5 of which were later rejected. In addition, the data 
on a large number of galaxies, present in the previous 
version of the sample, were re-measured. The majority 
of such data featured only slight modifications with re- 
spect to the previous version. This suggests that these 
data are of good quality. Nevertheless, we changed the 



data for 190 out of 1623 galaxies from the previous sam- 
ple. For 88 galaxies the changes were minor and did not 
drastically affect the results. However, the data for 47 
galaxies were changed significantly. This means that for 
these galaxies either the H i line width changed by more 
than 20kms~ or the radial velocity changed by more 
than 40kms^^, the prior being more common. The 
maximal change of velocity exceeded 2000 km s~^ and 
the maximal change of Hi line width was 126kms~ . 
We will give the details of the sample construction in 
Section ??. 

Thus, we can say that our sample has 37 added -t- 47 
replaced = 84 essentially different data, which equals to 
5% of the sample volume. This increase of the volume 
is smaller than in the previous versions of the sample. 
This is likely due to the fact that most galaxies in the 
field of regard of the Arecibo radio telescope were al- 
ready measured and the observations in the southern 
hemisphere are lagging behind. The progress in the 
Arecibo zone can be provided by the remaining 60% of 
the ALFALFA survey and the improved measurements 
of the rejected galaxies. The greatest impact would be 
provided by the southern sky observations, since these 
are quite scarce and the addition of new data would not 
only increase the sample volume but also make it more 
uniform and symmetric. 



3 Models of the collective motion 



In our previous articles ( IParnovskv et al.l l2001l: Par 
novsky & Tugay i2004 : .Parnovskv fc Parnowskil 



2010a 



20111 ) we described the models of the collective motion 
we used. Here we will only give two models given there. 
We will start from the most complex of the previously 
considered models, namely the DQO-niodel. 



V ^ R + V^'P + F'^'"" + V""* + AV. 



(1) 



Here ^ is a radial velocity of the galaxy in the CMB 
isotropy frame, R = Hr is the Hubble velocity, r is the 
distance to the galaxy, AV^ is a random error, V'^'^p, 
yqua yoct ^^^^ dipolc (D), quadrupolc (Q) and 
octupole (O) components of the non-Hubble cosmic 
flow. They are given by the following expressions: 

= RQ^kn,nk 
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Here we use Einstein's convention - summation by re- 
peated indices. Ui are the Cartesian components of the 
unit vector n towards the galaxy, connected with the 
galactic coordinates I and b in the following way: 



ni = = sin 6, 

= ~ cos I cos 6, 
n3 — Uy — sin /cos 6. 



(3) 



The dipole component (bulk motion) is described by the 
vector D. The quadrupole component (cosmic shear) is 
described by the symmetrical traceless tensor Q. It has 
5 independent parameters qi. The octupole component 
can be described by one rank 3 tensor, but we divide 
it into a trace characterised by vector P and a tensor 
O, which is antisymmetrical with respect to each pair 
of indices. The latter has 7 independent parameters o^. 

Hubble velocity is determined from the generalised 
Tully-Fisher relation in the "angular diameter -Hi line 
width" version by the following formula 



R^{C\+ C2B + C3BT + dU)^ 



(4) 



where W is the corrected H i line width in km s^ mea- 
sured at 50% of the maximum, a is the corrected ma- 
jor galaxies' angular diameter in arcminutes on red 
POSS and ESO/SERC reproductions, U is the ratio 
of major galaxies' angular diameters on red and blue 
reproductions, T is the morphological type indicator 
[T = It — 5.35, where /( is Hubble type; /t = 5 cor- 
responds to type Sc), and B is the surface brightness 
indicator {B ^ Isb—'^, where Isb is the surface bright- 
ness index from RFGC; brightness decreases from I to 
IV) . Note that the statistical significance of each term 
in eg. (Ill) is greater than 9 9% according to the F-test 



( Fisher 



19541 : lHudsonlll964r ). 

the DQO-model contains 



Thus, the DQO-model contains 24 parameters, 
namely 3 components of the vector D, 6 coefficients 
Ci, 5 parameters qi of the tensor Q, 3 components pi of 
the vector P and 7 parameters Oi of the tensor O. By 
rejecting V°''* we get a simpler DQ-model with 14 com- 
ponents. Further rejecting leads to the simplest 
D-model with 9 components. 



In the article (iParnovskv fc Parnowski 



base of the results of Kudrva fc Alexandrov 



2011) on the 



(|2002l . l2004l ) 



we also introduced relativistic models of motion based 
on the idea that for the homogeneous isotropic cosmo- 
logical models the dependence of the velocity V = cz 
on the angular diameter distance R = Hr in the next 
order in R has the form 



The coefficient 7 is expressed through the deceleration 
parameter q by 
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3 + g 
2c 



For the standard ACDM cosmology we have 
g = ^ - = -0.61, 



(6) 



(7) 



where flm and f^A are the relative densities of matter, 
including dark matter, and dark energy respectively. 
Numerical estimations are based on the resu lts of 7-year 
WMAP observations dKomatsu et al.lboiof ). Introduc- 
ing (O into ([U we obtain 
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3.98- 10~Skm"\ 



(8) 



In this article we use this so-called semi-relativistic 
model of galaxy motion with fixed 7 in the inhomoge- 
neous space-time, which has the form 



V° 



- 70-R^ + Ay. 



(9) 



The use of 70 mitigates the negative impact of mea- 
surement errors in the Hi line widths and the angular 
diameters. Please refer to the articles f Parnovskv & 



2010bl . I2OIII ) for details. Since we consider 



Parnowski 

the terms proportional to R^ separately, we should re- 
move the terms quadratic in distance from the gener- 
alised Tully-Fisher relation (jlj: 



R={Ci+ C2B + C3BT + CiU) 



W 



1 



(10) 



Note that all the remaining terms in this equation are 
inverse proportional to the angular diameter a. 

It is possible to add the next two multipoles to the 
DQO-model by constructing the DQOX-model 



V = R + V^'P + y"" + V""* + V''^"'-' + AV 



(11) 



and the DQOXT-model 

V ^ R + + y?"" + V°''^ + V'^"^ + V^^'^'> + AV. (12) 



V = R + -fR^. 



(5) 
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The velocity components corresponding to the 16-pole 
and the 32-pole have the form 



+90" 



(13) 



yhex ^ R^XijkiniTijnkni = {xin\ + X2nj 
+ X3n| + XAn\n2 + x^n\nz + xenfni 
+ X7n|n3 + x^n^ni + x<jrv^n2 + xiQn\n\ 
+ Xiinlnl + Xi2'n'ln\ + Xi^n\n2'n^ 
+ xunlniris + Xizn\nin2) , 

T/(32) ^ R'^Tijkim'riinjnkninrn = i?** (iinf 
+ t2fi2 + ^3^3 + t^n\n2 + t^nfri^ 
+ tQU^n^ + t^n^ni + tg,n\ni + tgn'^n2 
+ tionlnl + tiinlnl + tunlnl 
+ tianln^ + tunlnl + tiznln2 
+ tienln2n3 + tnnln^ni + tisnlnin2 
+ tign^rij^s + t2onlnlni + t2in\n\n2) ■ 



However, it is more correct to introduce these mul- 
tipoles to the semi-relativistic model of motion 

Vrei = R + V^'P + V'''^'' + V'"'^ + V^'''-' +-1oR^ + IW (14) 



and 



(15) 



where R is given by the relativistic TuUy-Fisher re- 
lation (ITU|) . This is due to the C'5 (■^) term in the 
non-relativistic Tully-Fisher relation ([4]), which causes 
different multipoles to be indistinguishable. 




180° 



Fig. 1 A fragment of the Mollweide projection of the ce- 
lestial sphere. Solid lines show the la, 2a and 3(7 confidence 
boundaries of the bulk flo w apices at fOOfc~^ Mpc. Black - 
the new sa mple, green - llParnovsky fc Parnowskill20f0a[) . 
sky blue - l|Parnovsky fc Tugav 2004h . The results of other 
authors are marked with red lo wercase lett ers, a - f Lvnden- 
BeU et al. 1988). b - llLauer fc Postman f 9941). c - (Hudso n 



et a l. Il995l). d - toale et al.lll999l). e - llDekel et al.llf999l) , 
f - llda Costa et aLlbOOOl ) g - (iParnovskv et al.ll200f[ ). h - 
llKudrva et al.ll2003l). i - (iHudson et al ll2004D. i - f Watkins 
et al. I2OO9I ), k - jP eldman et al.ll2010l) . 1 - the composite 
sample from (|Ma et al. 20f f ) . The apices in different models 
of motion are shown with uppercase letters. D - D-model, 
Q - DQ-model, O - DQO-model, X - DQOX model, T - 
DQOXT-model. For DQO-, DQOX- and DQOXT-models 
the non-relativistic versions are shown with black colour and 
the semi-relativistic - with blue colour. The numbers show 
the directions of the quadrupole eigenvectors in the DQO- 
model. f - maximum, 2 - intermediate, 3 - minimum. 



4 Results 

We process the new sample taking into account each 
model of motion described above. We created sub- 
samples with a limitation on the distance in the non- 
relativistic D- model R < Rmax- This is the most 
straightforward way for limiting the sample depth. The 
mean distances for the two most often used subsam- 
ples with Rmax = 80h~^ Mpc and Rmax — I00h~^ Mpc 
are equal to 45Ah~^ Mpc and 51.5h^^ Mpc respectively. 
For each subsample we determined the regression coef- 
ficients using the least square method simultaneously 
in the Tully-Fisher relation and the model of motion. 
The results obtained allow estimating the quality of the 
sample. 

The main change since the previous sample is the 
slightly different apex of the dipolar component. This 
is demonstrated in Figure[T]where the boundaries of la, 
2(7 and 3(7 confidence areas of the apices of the dipolar 
component are shown for 3 different samples using the 
D-niodel with Rmax = 100ft.~^Mpc. The new sample 
is shown with black lines, the sample introduced in the 



articles ([Parnovskv fc Parnowskill2Q09l . l2010a| ) is shown 
with light green lines, an d the sample introduced by 
Parnovskv fc Tugavl (|2004r ) is shown with sky blue lines. 
The Figure [1] is a Mollweide projection of a part of the 
celestial sphere. 

ft can be clearly seen that the apex travels in the 
same direction with the improvement of the sample. 
As a result, the apex deviated from the results of most 
authors (red lowercase lett ers in Figure [J) and be - 
came closer to the results of iLauer fc PostmanI (jl994l ). 
However, the obtained value of the bulk flow velocity 
278kms~^ is much smaller then that given bv Lauer 
fc Postman (Il994 ) and is consistent with ACDM cos- 
mology. Naturally, we analyzed the reasons behind this 
phenomenon. This change of the apex is mostly due to 
the improved measurements of the H i line widths of the 
galaxies, which were already present in the sample. All 
such re-measurements were double-checked by us and 
appeared to better match the regression relation. It is 
important to note that the changes in the H i line width 
have approximately equal probabilities to be positive or 
negative, thus imposing no systematic shift. 
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Now let us discuss the sample compilation. This step 
is very important, because it is one of the few manual 
operations in the routine. Errors at this step can lead 
to underestimation of the quadrupole and higher mul- 
tipoles due to selection. 

The new data on radial velocities V and Hi line 
widths W were compared to the data from the pre- 
vious versio n of the sample described by Parnovskv 
& Parnowski (l2010af) and published in ( Parnovskv & 
Parnowskii20^~ The data about newly added galax- 
ies were included automatically. The data, whose differ- 
ence from the old ones either in or in was less than 
a few percent, were considered to be the same as before 
and were not changed. The data, whose difference from 
the old ones exceeded few percent, were considered as 
alternatives to the old data. Thus, we have a problem: 
if we have several alternative sets of data for a galaxy, 
which one we should choose? 

We addressed this problem in the following way. In 
Figure [2] we plotted all the galaxies with known H i 
line widths. On the horisontal axis we plotted the dis- 
tance according to D-model. On the vertical axis we 
plotted the difference between the redshift velocity and 
the value of the velocity following from the D-model 
at Rmax — lOOft."^ Mpc. The galaxies which enter the 
sample are marked with dots, and the rejected galaxies 
- with circles. Among all possible alternatives we chose 
those with smallest deviations. However, this cannot be 
done in a straightforward way because the value of the 
deviation depends also on the regression parameters, 
which are sample-specific. 

For this reason we calculated the deviations using 3 
different regression relations. The first one is based on 
the old sam ple fits parameters are given in (^ Parnovskv 
& Parnowski |20 1 Oaf ) ) . The second one is based only on 
the galaxies, which enter the regression and have only 
one set of data. The third one initially includes all the 
candidate data, i.e. all the different sets of data for 
each galaxy were included as separate entries. Then 
we start rejecting the galaxies from the most obvious 
outliers. At each stage the regression relation is recal- 
culated. However, we must stop at some point. Note 
that we included in the sample some galaxies with the 
deviations up to 3.5cr. Could we reject more galaxies 
and reduce the standard deviation of our sample? Let 
us demonstrate that this will lead to selection. 

There are two effects which cause the deviation from 
the D-model. The first one is random scatter, which we 
will roughly estimate with a normal distribution with 
zero mean and variance. The second one is the im- 
pact of the quadrupole and higher multipoles, which we 
will denote as A. Note that A greatly varies across the 
sample, taking the largest values at the edges of sample 
and in the vicinity of attractors. 



Which galaxies will enter the sample and which will 
be rejected? If we assume the deviation from the D- 
model to be purely random, i.e. A = 0, we will reject 
the galaxies whose deviations exceed some threshold. 
For the 2a threshold the percentage of rejected galaxies 
is 2.3% for positive deviations and the same for negative 
ones. For the 2.5cr threshold these percentages are 0.6% 
in both ways, and for the 3(T threshold they are both 
less than 0.1%. Rejection of these galaxies yields no 
systematics. 

Now let us consider a region of space where the 
higher multipoles account for a constant impact of 
|A| = a. The percentage of galaxies in such areas can 
be roughly estimated as 10% of the sample. A non-zero 
value of A will lead to an effective shift of the thresholds 
by A, positive for the deviations with the same sign as 
A and negative in the opposite case. This effect causes 
a selection which leads to an underestimation of higher 
multipoles. In this case the percentages of the rejected 
galaxies will become asymmetric. For the 2a threshold 
the percentages of rejected galaxies are 15.8% in the 
direction of A and less than 0.1% in the opposite one. 
For the 2.5a threshold the percentage ion the direction 
of A is 6.5%, and for the 3a threshold it becomes 2.3%. 
The percentages of rejects in the opposite direction for 
2.5cr and 3a thresholds are negligible. Thus, the intro- 
duced selection for this particular area is 16% for the 
2a threshold, 6.5% for the 2.5cr threshold, and 2.3% for 
the 3(7 threshold. Multiplying this by 10% of the sam- 
ple volume (about 170 galaxies), we get the expected 
numbers of erroneously rejected galaxies: 27 galaxies 
for the 2a threshold, 11 galaxies for the 2.5cr threshold 
and only 4 galaxies for the 3a threshold. 

To verify these simple estimations we considered the 
subsamples with different deviation thresholds. They 
are shown in Figure [2] with horizontal lines correspond- 
ing to 2000, 2500 and 3000 kms"^ thresholds. Combin- 
ing these thresholds with DQ- and DQO-models and 
4 different limits on distance, namely 80h~^, 100h~^, 
lAOh^^ and 170/i^^ Mpc, we get 32 different combi- 
nations. For each of them we calculated the velocity 
of the dipolar component as well as the maximal Ai 
and the minimal A3 eigenvalues of the shear tensor Q 
(the intermediate eigenvalue is not independent and is 
equal to A2 = — Ai — A3). The maximal dipolar ve- 
locity D among all possible combinations was equal to 
304kms~ , which is consistent with the ACDM cos- 
mology. The values of the dipolar velocity and the 
eigenvalues are given in Table [1] One can see that 
the quadrupole significantly drops with the 2500 kms~^ 
threshold for all subsamples. For the DQ-model with 
Rmax = 80h~^ Mpc subsample it drops already at the 
3000 kms~^ threshold. The selection criterion we used 
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Fig. 2 The deviation from tlie D-model vs. tlie radial 
distance for all galaxies. Horizontal lines show the 2000, 
2500 and 3000 kms~^ thresholds 



to construct the sample appeared to be free from such 
problems. 

We obtained the norm of the dipolar velocity compo- 
nent for different models and different subsamples. In 
Figure[3]we plotted its dependence on the sample depth 
Rmax in the framework of D, DQ and DQO models. In 
the simplest D-model the norm grows up to the distance 
100/i~^ Mpc, in the DQ-model it is almost constant be- 
yond 70h~^ Mpc, and in the DQO-model it has a mini- 
mum at about 82.5h~^ Mpc (the minimal norm is non- 
zero at 1.5(7 confidence level). For Rmax > 120/i^^ Mpc 
these values are almost constant due to the small num- 
ber of more distant galaxies in our sample (47 out of 
1661). 

These results are obtained for unit weights of all 
data points. Let us see how these results will be af- 
fected by weighting. Let us choose the weights, which 
decrease with distance inverse proportional to A^{R), 
where A(i?) is a maximum deviation of the velocity 
from the D-model at the distance R. After calculating 
the norm of the dipolar velocity component for the sub- 
samples limited at 80, 100 and 120/1"^ Mpc we obtain 
for the D-model 243, 251, and 243knis^^ respectively, 
for the DQ-model 245, 248, and 244kms~\ and for the 
DQO-model 250, 262, and 280 kms"^ Ah these values 
do not differ essentially from those with unit weighting 
and still are consistent with ACDM cosmology. 

Additionally, we analyzed the impact of the 16-pole 
and 32-pole on the velocity field. The goal of this anal- 
ysis was to analyze how the inclusion of higher mul- 
tipoles aff'ects the dipole, the quadrupole and the oc- 
tupole. The statistical significance of both multipoles 




8000 
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Fig. 3 The dependence of the norm of the dipolar velocity 
component on the sample depth Rmax in the framework of 
D, DQ and DQO models 



exceeds 99.9% according to F-test and the inclusion of 
the 16-pole reduces the RMS error by 4.5%. The semi- 
relativistic DQOX-model gives a much better apex as 
compared to the semi-relativistic DQO-model and the 
value of the dipolar velocity 236 km s^ is also reason- 
able. Thus, the semi-relativistic DQOX-model can be 
used for deep samples. The DQOXT-model contains 
too much parameters and it is too early to be considered 
unless the sample volume and accuracy is essentially im- 
proved. The maximal dipolar velocity in these models 
appeared to be equal to 342kms~^ in the DQOXT- 
model with Rmax — SOh~^ Mpc, which is marginally 
consistent with ACDM cosmology. 

In FigurelUwe plotted the apices of the dipolar veloc- 
ity in different models of motion. The letter D denotes 
the D-model, Q - the DQ-model, O - the DQO-model, 
X - the DQOX model and T - the DQOXT-model. 
For DQO-, DQOX- and DQOXT-models both the non- 
relativistic and the semi-relativistic versions are shown 
with black and blue colour respectively. One can see 
that the non-relativistic DQO-model gives the apex 
closest to the results of most authors. This is consis- 



tent wilh_^^ earlier results (jParnovskv fc Parnowski 
20081 . 2011 ) of Monte Carlo simulations, which demon- 
strate that the DQO-model gives the best fit to the 
data. However, the apex in the semi-relativistic DQO- 
model is way beyond all reasonable confidence bound- 
aries. This is due to the low value of the dipolar velocity 
in this model, which yields large errors in the determi- 
nation of the apex. 
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Table 1 The impact of the deviation threshold on the dipolar velocity and the maximal and minimal eigenvalues of the 
quadrupole 



Max deviation, km s ^ 


D, kms ^ 




Ai, 


% 




A3, % 


Full 


3000 


2500 


2000 


Full 


3000 


2500 


2000 


Full 


3000 


2500 


2000 


Rmax, /i^^Mpc 


DQ-modcl 


80 


217 


197 


186 


195 


7.4 


3.1 


2.5 


3.5 


-5.8 


-3.7 


-4.0 


-3.3 


100 


224 


213 


195 


227 


6.2 


6.2 


3.5 


2.5 


-4.2 


-4.0 


-3.3 


-3.1 


140 


233 


222 


196 


225 


3.7 


3.8 


2.4 


2.8 


-3.9 


-4.5 


-3.4 


-4.0 


170 


228 


225 


196 


227 


4.0 


3.9 


2.4 


2.9 


-4.4 


-4.5 


-3.4 


-4.2 


Rmax, /l~^MpC 


DQO-model 


80 


158 


304 


269 


258 


7.1 


7.0 


5.6 


3.5 


-5.6 


-6.5 


-6.3 


-4.0 


100 


178 


212 


283 


275 


7.3 


7.1 


4.7 


2.9 


-4.8 


-4.8 


-4.4 


-4.0 


140 


292 


261 


269 


236 


7.0 


7.3 


5.2 


4.3 


-6.5 


-6.5 


-5.5 


-6.0 


170 


304 


260 


269 


237 


7.1 


7.3 


5.2 


4.3 


-6.7 


-6.4 


-5.5 


-6.2 



For reference we give the parameters of the non- 
relativistic DQO-model at Rmax — lOO/i^^Mpc. The 
residual mean square error is a = 1112 km s^"'^. 

The coefhcients of the generalised Tully-Fisher rela- 
tion are 



Ci = 18.0 ± 1.3, Ca = 1.59 ± 1.8, 

C3 = -0.27 ± 0.11, C4 = 6.5 ± 1.1, 

C5 = (-7.2 ± 1.1) • 10-^ Ce = -919 ± 82. 



(16) 



The dipolar velocity is equal to D ^ 178 km s^^ and 
points to the apex with the galactic coordinates I — 
303°, 6 = -1-11°. The components of the dipolar velocity 



are equal to — 35 ± 74 kms , D. 
D„ = -147 ± 94 kms" 



Feldman et al 



94 ±95 k ms"', 
(|2010f ) using 



a model similar to DQO on the sample with the same 
depth 100h~^ Mpc obtained the dipolar velocity D = 
416 ± 78kms"^ pointed towards / 282 ± 11°, b = 
+6 ±6°. 

The quadrupolar component is described by its 5 
irreducible components: 



qi = (7.2 ±1.5)%, 92 - 
q3 = (-0.8 ±2.0)%, 94 
95 = (-1.6 ±2.6)%. 



(-2.7 ±1.6)%, 
= (1.9 ±2.4)%, 



(17) 



It can be also represented by its eigenvalues Ai = 
(7.27 ± 1.54)%, A2 = (-2.43 ± 1.46)%, A3 = (-4.84 ± 
1.47)% and eigenvectors pointing respectively to / = 
118°, b = +85° (Canes Venatici); / = 341°, b = 
+4° (Scorpius/Norma); and I = 71° , b = -4° 
(Cygnus/Volans) and in the opposite directions: Sculp- 
tor, Auriga and Centaurus/Vela/Carina. The direc- 
tion of the eigenvector corresponding to the maximum 
eigenvalue is close to the Supergalactic plane and the 
direction corresponding to the intermediate eigenvalue 
is not far from the direction towards the Great At- 
tractor. These directions are marked in Figure [T] with 



numbers 1 (maximum), 2 (intermediate) and 3 (mini- 
mum). It is inter esting to compa r e the se values with 
those obtained by Courtois et al.l ( 2012h . who consid- 
ered a sample of 1797 galaxies within 30h~^ Mpc. Their 
directions of the eigenvectors demonstrated good agree- 
ment with our results for 100/i~^Mpc. However, our 
results for 30h~^ Mpc were in not so good agreement 
due to a small number of galaxies in this subsample. 

The octupolar component is described by 10 irre- 
ducible components: 

Pi = (2.2 ±2.1) • 10-6skm"\ 
P2 = (1.0 ±2.5) • 10-Skm"\ 
P3 = (0.4 ±3.1) • 10-6skm"\ 

01 (3.0 ± 1.5) • 10-6skm~\ 

02 = (1.0 ± 1.8) ■ 10-'^skm-\ 

03 = (6.7 ± 1.8) • 10-6skm-\ 

04 ( 



(-2.9 ±2.2) • 10-Skm"\ 

05 = (5.5 ±2.0) • 10-6skm"\ 

06 = (-2.0 ±2.1) ■ 10'Skni"\ 

07 = (22.0 ±7.3) • 10-6skm~\ 



(18) 



The trace vector P, which is a part of the reduced oc- 
tupole tensor and sho ws how the dipolar com ponent 
changes with distance ( Parnovskv et al]|200lh . points 
towards the direction I = 23°, b = +64° (Bootes). 

In Figure |4] we plotted the radial component of the 
velocity field calculated in the non-relativistic DQO- 
model with above parameters at i? = 80h~^ Mpc. The 
global minimum lies in the direction opposite to the 
apex and the global maximum matches the direction 
of the maximum eigenvector of the quadrupole com- 
ponent. This can explain the shift of the apex in the 
D-model towards the latter. 

Note that each multipole of order n > 2, i.e. quadru- 
pole and higher, can be represented in a reduced form, 
in which the multipoles of orders n — 2, n — 4 and so 
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-90° 



Fig. 4 The radial component of tlie velocity field calcu- 
lated in the non-relativistic DQO-model at J? = 80h~^ Mpc 



to be statistically significant. The inclusion of the 16- 
pole was beneficial for the semi-relativistic model and 
improved the obtained results. 
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on are separated. However, it is important to realize 
that these lower order multipoles have the same depen- 
dence on the radial distance as the initial one. We use 
the reduced form for the quadrupole, separating the 
Hubble constant H, and for the octupole, separating 
the trace vector P, which describes t he change of the 



bulk flow velocity with distance (see (jParnovskv et al. 



20011 ) for details). We do not use the reduced form for 



the 16-pole and the 32-pole, because their separated 
lower-order parts do not have obvious physical mean- 
ing. In addition, since these lower-order parts have a 
different dependence on the distance than the quadru- 
pole and the octupole, there is no interference between 
them and both their reduced and full forms can be used 
interchangeably. 



5 Conclusion 



We described a new sample of RFGC galaxies intended 
for the study of large-scale cosmic flows. A preliminary 
verification with the methods described in our previous 
articles has shown that this new sample can be used as 
a basis for the study of collective motions of galaxies. 
We plan to modify our routine to reconstruct the dis- 
tribution of matter density simultaneously with the pa- 
rameters of the model of collective velocity field. In this 
article we applied the old routine to the new sample to 
verify its quality. The obtained results show principal 
agreement with the previous ones. It is essential that 
in all considered cases we obtained the values of the 
dipolar velocity, which are consistent with the ACDM 
cosmology. Once again, it was demonstrated that the 
samples with significant depth can not be considered in 
the framework of a simple D-model featuring only Hub- 
ble expansion and the bulk flow. For our sample with 
the depth 80 — 100/i^^Mpc the optimal choice is the 
DQO-model, which also takes into account the cosmic 
shear and the octupolar components. For the first time, 
higher-order multipoles were considered and appeared 



Yet another sample of RFGC galaxies 
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